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We present an exploratory calculation of the standard model AS — 2 matrix element relevant for indirect CP 
violation in A' ^ tttv decays. The computation is performed with overlap fermions in the quenched approximation 
at /3 = 6.0 on a 16'' x 32 lattice. The resulting bare matrix element is renormalized non-perturbatively. Our 
preliminary result is B^°^{2 GeV) = 0.61(7), where the error does not yet include an estimate of systematic 
uncertainties. 



1. Introduction 

j{0 _ ^0 niixing induces indirect CP violation 
in K ^ TTTT decays and is governed, in the stan- 
dard model, by the matrix element 
16 
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where C'as'=2 is the four-quark operator given be- 
low (Eq. (2)) and where Bk parametrizes devia- 
tion from the vacuum saturation approximation. 

The benefit of using overlap fermions [1] 
to compute this matrix element is their exact 
SU{Nf)L X SU{N-f)R flavor-chiral symmetry at 
finite lattice spacing [2]. This implies that the 
AS* = 2 operator is simple to construct on the 
lattice (unlike with staggered fermions) and does 
not mix with chirally dominant operators [3] (as 
it does with Wilson fermions). Overlap fermions 
further guarantee full O (a)-improvement. On a 
different note, this exploratory computation tests 
the feasibility of calculating matrix elements of 
four-quark operators with overlap fermions. 

2. Simulation details 

We use Neuberger's overlap action [1] 
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with 



D = + X/VJOX^ and X = Dw - p/a , 

where Dw is the massless Wilson-Dirac operator. 

The renormalized AS = 2 operator of Eq. (1) 
is related to the bare lattice operator through 



(2) 



where the lattice operator with correct chiral 
properties is 

OTsU = [s7m(1 - 75)d][s7M(l - 75)d] , 

withd= {l-aD/2p)d. 

The computation is performed with the Wil- 
son gauge action in the quenched approximation 
at /3 = 6.0 (a-i(ro) ~ 2.1 GeV) on a 16^ x 32 
lattice with a statistics of 80 configurations and 
with p = 1.4. The pseudoscalar mesons simulated 
are composed of degenerate quarks with masses 
ruq = 0.04, 0.055, 0.07, 0.085, 0.1 corresponding to 
a range from ~ m<;/2 — > 1.3ms. Details regard- 
ing the implementation and the inversion of the 
overlap operator can be found in [4] . 

The three-point function used to extract the 
AS" = 2 matrix element and its asymptotic time 
behavior are 
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with J^=d7^(l 



We allow for a time- 



reversed contribution because our lattice is rather 
short in the time direction and this contribution 
can be as large as 3-4% for our lightest quark 
mass at = T — j/o = 8a, a point which is 
included in our fits. In Eq. (3) S measures the 
finite-volume shift in the energy of the two kaon 
state. This shift varies from about 2% to 5% in 
going from our heaviest to our lightest mass and 
is consistent with zero within roughly three to one 
standard deviations. Rk is clearly intimately re- 
lated to {0\O''^^'L2\KK), but finite volume effects 
can be significant and we choose not to pursue 
the study of this time- reversed matrix element. 
The fit that we retain for determining the for- 
ward matrix element is to times xq/cl and yo/a in 
the intervals 6^8 and 24 26 with yo > xq. To 
obtain K^j J,^|0)| and Mk we compute two-point 
functions of left-handed currents which we sym- 
metrize in time and fit in the time range 5 — > 15. 

The motivation for using left-handed currents 
is that they eliminate spurious zero-mode contri- 
butions which appear in finite volume and in the 
quenched approximation because the overlap op- 
erator satisfies an index theorem [5], as a good 
Dirac operator should. In our correlators, 
topological zero modes, 4>o{x), contribute as 

(j)o{x)(f>l{y) 

7m 1 - 75) 7m 1 - 75) , 

ruq 

which vanishes because of the chirality of these 
modes. The use of left-handed currents is pro- 
posed in [6] where it is suggested that the weak 
interactions of Goldstone bosons be studied in the 
e-regime of Gasser and Leutwyler. 

3. Chiral behavior of Bk 
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where Jq is the conjugate of . Quenched chiral 
perturbation theory (QxPT) predicts [7] 
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Figure 1. Chiral behavior of Bk- Fits are to Eq. 
(4), as described in the text. The vertical dashed 
line indicates the physical point. 

where B, F and M are the leading order val- 
ues of the i3-parameter, leptonic decay constant 
and pseudoscalar mass, respectively. is a scale 
which corresponds to O (p*) contact terms in the 
effective theory and the term proportional to b is 
added to parametrize likely higher order contribu- 
tions. We fit this functional form to our results for 
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6 = or not. The fit parameters are B , Ab 
and b when it is not fixed to zero. ^ The results 
of these fits are shown in Fig. I. We find that the 
mass-dependence of Bk is consistent with expec- 
tations from xPT , though results at lower masses 
and with better statistics would be necessary to 
confirm the coefficient of the logarithm. Until 
such results become available, we choose not to 
quote a number in the chiral limit. At the physi- 
cal point Mk = M^^^ , however, where the value 
of the i3-parameter obtained is insensitive to the 
choice of functional form, these considerations are 
not important. 

4. Non-perturbative renormalization 

We perform all renormalizations non- 
perturbatively in the RI/MOM scheme a la [8]. 
Thus we fix gluon configurations to Landau gauge 
and compute numerically approriate, amputated 
forward quark Green's functions with legs of mo- 



mentum p : 



Then we use the fact that with 



^Note that the definition of b changes with F. 

''The result for Bk given below is obtained from the chiral 

fit with F = FP''^" and fe ^ 0. 
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Figure 2. Scale dependence of the renormalization 
constants for Bk defined in the text. The three 
curves for Z^^'^'^{p) are obtained by allowing the 
scale parameter in the coupling of [9] to take the 
values A^J2, A^^ and 2Ap^. 

overlap fermions the renormalization constants 
ZjL = Zv and Zas=2 = Zvv+aa to define 
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where To is the value of the non-perturbative, 
amputated Green's function of operator O pro- 
jected onto the spin-color structure of O. 

We find that in a range from approximatively 
1.5 to 2 GeV the p-dependence of Zg^^ matches 
2-loop running obtained by combining the results 
of [10] and [11], with Nf = and as taken from 
[9]. This is shown in Fig. 2 where we plot, vs p, 
ZW'^ which we define as the ratio of ZW to the 
2-loop running expression. Also shown in Fig. 2 
are the non-perturbative Zg^^ (p) and its 1-loop 
counterpart obtained by expanding in 

Us the ratio of ^^5=2^) i^A^'^^)^- The lat- 
ter are taken from [12] and are properly matched. 

The value of ZS*^^ that we use is the result of fit- 
is a- 

ting our data to a constant, in the momentum 
range shown in the plot. The renormalization 
constant in the MS — NDR scheme is then com- 
puted through multiplying this value by the 2- 
loop running expression of [10] with Nf and 
chosen as above. 

5. Conclusion 

Though costly, it is entirely feasible to com- 
pute weak matrix elements with overlap fermions. 



at least in the quenched approximation. The ad- 
vantages of this approach are the chiral symmetry 
and the O (a)-improvement. 

Our preliminary results are Z^'^^ ~ 1.29(1) 
for the renormalization constant and, for the B- 
parameter: 



rRGI 



0.85(10) 



rNDR 



(2 GeV) = 0.61(7), 



where the errors do not yet include an estimate of 
systematic uncertainties. For the S-parameter, 
agreement with continuum-limit world averages 
(see e.g. [13]) based on the staggered result of [14] 
is excellent. Agreement with the results of [15], 
obtained with a different formulation of overlap 
fermions, perturbative matching and larger quark 
masses is good. Finally, our result is respectively 
a half and one standard deviations above the re- 
sults of [16] and [17], both obtained using domain- 
wall fermions with a finite fifth dimension. 
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